We used the fPrioritizrRun() function to streamline the runs. fPrioritizrRun() arranges all the data and makes sure that they’re in the right format, then runs problem() and solve(). This function also saves the summaries each plan (cost, % area protected and number of representative features that achieved 30% or more protection).
The relevant scripts for this step is found in ~/scripts/ with the prefix 10_Prioritizr.
# the function is in:
source("scrpits/10a_Prioritizr_fPrioritizr.R")
# to run the function with climate-smart data
source("scripts/10h_Prioritizr_fPrioritizrIUCNSmart_NoProv.R")
# to run the function with climate-uninformed data
source("scripts/10i_Prioritizr_fPrioritizrIUCNUninformed_NoProv.R")
# to create no-regret closures with the climate-smart plans
To create no-regret plans, we used thre function fCreateNoRegret() which intersected the selected PUs in all three climate scenarios. The relevant scripts for this step is found in ~/scripts with the prefix 11_NoRegret.
# the function is in:
source("scripts/11_NoRegret_fCreateNoRegret.R")
# the runs are found in:
source("scripts/11e_NoRegret_fCreateNoRegretIUCNRuns_NoProv.R")
As mentioned in the previous section, we will only be showing results of Runs 2 (Max = 0.9, Min = 0.1), 5 (Max = 0.6, Min = 0.1) and 8 (Max = 0.3, Min = 0.1).
SSP1-2.6 run
SSP126_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_SSP126.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP2-4.5 run
SSP245_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_SSP245.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP5-8.5 run
SSP585_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_SSP585.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
noregret_target90 <- readRDS("outputs/11_NoRegret/11e_IUCNRuns_NoProv/noregretclosures_Target90.rds")
uninformed_target90 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/02_Target90/solution_uninformed.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 1 x 12
## target_scenarios sum_area percent_area num_pu total_cost median_velocity
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Target90_summaryuninf… 6.72e7 79.0 25207 767309. 67242193.
## # … with 6 more variables: median_tvelocity <dbl>, median_RCE <dbl>,
## # median_RCEtvalue <dbl>, represented_features <dbl>, target <dbl>,
## # scenario <chr>
SSP1-2.6 run
SSP126_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_SSP126.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP2-4.5 run
SSP245_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_SSP245.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP5-8.5 run
SSP585_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_SSP585.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
noregret_target60 <- readRDS("outputs/11_NoRegret/11e_IUCNRuns_NoProv/noregretclosures_Target60.rds")
uninformed_target60 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/05_Target60/solution_uninformed.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 1 x 12
## target_scenarios sum_area percent_area num_pu total_cost median_velocity
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Target60_summaryuninf… 4.48e7 52.7 16807 393372. 44834353.
## # … with 6 more variables: median_tvelocity <dbl>, median_RCE <dbl>,
## # median_RCEtvalue <dbl>, represented_features <dbl>, target <dbl>,
## # scenario <chr>
SSP1-2.6 run
SSP126_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_SSP126.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP2-4.5 run
SSP245_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_SSP245.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
SSP5-8.5 run
SSP585_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_SSP585.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 0 x 12
## # … with 12 variables: target_scenarios <chr>, sum_area <dbl>,
## # percent_area <dbl>, num_pu <dbl>, total_cost <dbl>, median_velocity <dbl>,
## # median_tvelocity <dbl>, median_RCE <dbl>, median_RCEtvalue <dbl>,
## # represented_features <dbl>, target <dbl>, scenario <chr>
noregret_target60 <- readRDS("outputs/11_NoRegret/11e_IUCNRuns_NoProv/noregretclosures_Target60.rds")
uninformed_target30 <- readRDS("outputs/10_Prioritizr/10h-i_IUCNRuns_NoProv/08_Target30/solution_uninformed.rds")
cost is in USDperc_area_protected is the percentage of the Pacific ABNJ protectedprotected_PU is the number of features that have >= 30% (by area) protection## # A tibble: 1 x 12
## target_scenarios sum_area percent_area num_pu total_cost median_velocity
## <chr> <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Target30_summaryuninf… 2.24e7 26.3 8407 149748. 22426513.
## # … with 6 more variables: median_tvelocity <dbl>, median_RCE <dbl>,
## # median_RCEtvalue <dbl>, represented_features <dbl>, target <dbl>,
## # scenario <chr>